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Summary. We reconstruct the initial two-body relaxation time at the half mass 
radius for a sample of young 300 Myr star clusters in the large Magellanic cloud. 
We achieve this by simulating star clusters with 12288 to 131072 stars using direct 
iV-body integration. The equations of motion of all stars are calculated with high 
precision direct iV-body simulations which include the effects of the evolution of 
single stars and binaries. 

We find that the initial relaxation times of the sample of observed clusters in the 
large Magellanic cloud ranges from about 200 Myr to about 2 Gyr. The reconstructed 
initial half-mass relaxation times for these clusters has a much narrower distribution 
than the currently observed distribution, which ranges over more than two orders 
of magnitude. 

1 Introduction 

The early evolution of young clusters of stars is of considerable interest because 
of the rich observational data which became available recently. Studying young 
star clusters is also important in order to obtain a better understanding of 
the conditions under which clusters are born. The initial conditions of star 
clusters have been debated actively over the years, but no consensus has been 
reached either by observing or simulating the formation process of stellar 
conglomerates. The main parameters which characterize a star cluster at birth 
(and any time later) are the number of stars, the mass function and the 
concentration of the stellar distribution. 

Part of the problem in determining the clusters' initial conditions stems 
from our static view of the universe, our inability to run simulations backward 
with time and our lack of understanding of the physics of the star(cluster) 
formation process. In this paper we approach this problem by starting with 
a series of simulations in which we fix most of the initial conditions but relax 
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the initial relaxation time. We subsequently calculate the evolution of the 
star cluster with time. With this series of controlled numerical experiments 
we trace back the observed parameters for a number of young star clusters in 
the large Magellanic cloud along parallel trajectories in parameter space. 

2 Simulations 

We focus on young t ^ 300 Myr star clusters, because excellent observational 
data is available for a number of these P^. For the simulations we adopt 
the Starlab software environment which acquires it's greatest speed with 
GRAPE-6 special purpose hardware [3J0] 1 . Our simulations are performed 
on the GRAPE-6 hardware at the University of Tokyo and the MoDcStA 2 
platform at the University of Amsterdam. 

The simulated star clusters are initialized by selecting the number of stars, 
stellar mass function, the density profile, binary fraction and their orbital ele- 
ments. For our most concentrated model (simulation #1) we adopt the initial 
conditions derived by Portegies Zwart et al. [5] to mimic the 7-12 Myr old 
star cluster MGG-11 in the star-burst galaxy M82. In this paper, however, 
we extend the evolution of this simulated cluster to about 100 Myr &. Sub- 
sequent simulations are performed with a larger cluster radius, resulting in a 
longer initial relaxation time (for details on the simulation models #{1..4}c 
with 128k stars as listed in Tab^and see 0). The stellar evolution model 
adopted is based on |Jj , and the binaries are evolved with SeBa [Hj • 

We summarize the selection of the initial conditions for simulation $dc 
(see Tab.n see also 0): first we selected 131072 stars distributed in a King 
[§] density profile with Wo = 12 and with masses from a Salpeter initial 
mass function (x = —2.35) between 1 M Q and 100 M & . The total mass of the 
cluster is then M ~ 433000 M@. The location in the cluster where the stars 
are born is not correlated with the stellar mass, i.e. there is no primordial mass 
segregation. Ten percent of the stars were randomly selected and provided with 
a companion (secondary) star with a mass between 1 Mq and the mass of the 
selected (primary) star from a flat distribution. The binary parameters were 
selected as follows: first we chose a random binding energy between E — 10 kT 
(corresponding to a maximum separation of about 1000 Rq). The maximum 
binding energy was selected such that the distance at pericenter exceeded four 
times the radius of the primary star. At the same time we select an orbital 
eccentricity from the thermal distribution. If the distance between the stars 
at pericenter is smaller than the sum of the stellar radii we select a new semi- 
major axis and eccentricity. If necessary, we repeat this step until the binary 
remains detached. As a result, binaries with short orbital periods are generally 
less eccentric. We ignored an external tidal field of the Galaxy, but stars are 
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removed from the simulation if they are more than 60 initial half-mass radii 
(jhm) away from the density center of the cluster (100 ?\ m for the 12k models, 
see Tab-HJ. 

For the other simulations #2c, #3c and #4c, we adopt the same realization 
of the initial stellar masses, position and velocities (in virial N-body units |10|) 
but with a different size and time scaling to the stellar evolution, such that 
the two-body relaxation time (i r ix) for simulation #2c is four times that of 
#lc, for simulation #3c we used four times the two-body relaxation time of 
what was used for simulation #2c, etc. the initial conditions are summarized 
in Tab.H 

We subsequently generate additional initial realizations for clusters with a 
smaller number of stars (12k stars for models #{1..4}a versus 128k for models 
#{1..4}c). To study the effect of the initial density profile on the results are 
perform another set of simulations with 12k stars but with a Wq = 6 initial 
density profile rather than the highly concentrated Wo = 12. The simulations 
with N = 12k are constructed without primordial binaries. 

Table 1. Conditions for the performed calculations with a range of number of stars 
(lk = 1024 stars) and cluster virial radii. The columns give the model name, the 
number of stars in the initial model, the concentration parameter (Wo), the initial 
half-mass relaxation time (tl lx ), the initial virial radius and core radius and finally 
the initial crossing time. 



Run 


N 


Wo 


f l 

L rlx 




7*core 


tch 








[Myr] 


[PC] 


[PC] 


[Myr] 


#la 


12k 


12 


73 


2.27 


0.071 


0.258 


#lb 


12k 


6 


81 


2.77 


0.670 


0.352 


#lc 


128k 


12 


80 


1.27 


0.010 


0.032 


#2a 


12k 


12 


304 


5.72 


0.072 


1.05 


#2b 


12k 


6 


310 


6.98 


1.93 


1.010 


#2c 


128k 


12 


320 


3.20 


0.026 


0.129 


#3a 


12k 


12 


1200 


14.6 


0.151 


4.17 


#3b 


12k 


6 


1280 


17.8 


4.28 


4.19 


#3c 


128k 


12 


1300 


8.1 


0.066 


0.516 


#4a 


12k 


12 


4990 


36.3 


0.317 


16.6 


#4b 


12k 


6 


4920 


44.2 


14.9 


16.5 


#4c 


128k 


12 


5100 


20.0 


0.162 


2.07 
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After initialization we synchronously calculate the evolution of the stars 
and binaries, and solve the equations of motion for the stars in the cluster. 
The calculations are continued to an age of about 100 Myr. 

3 Results 

The main difference between simulations are the number of stars and since we 
ignore the external tidal field the number of stars drops only slightly during the 
simulations, but the half-mass relaxation time increasing substantially from 
its initial value. At the same time the cluster structure changes by becoming 
less concentrated. The latter is mainly attributed to mass segregation and 
stellar mass loss (see also ^2 El EH El El)- I n El we compared some of the 
characteristics of a subset of the here presented simulations with the observed 
sample of young star clusters in the LMC. Here, in this proceedings paper, we 
limit ourselves in comparing the currently observed two-body relaxation time 
and compare these with the results of our simulations. 

In figure^we present the evolution of the two-body relaxation time at the 
half-mass radius for the simulated clusters #l{a..c} to #4{a..c}. The relax- 
ation time is roughly constant for the first ~ 3 Myr for simulations #l{a..c}, 
=#=2{a..c}, #3c and #4c, to increase at later time. For simulations #3a and 
$=3b the half-mass relaxation time starts to increase only after about 6 Myr, 
and after about 16 Myr for simulation #4a and #4b. We attribute this ef- 
fect to the relatively slow response of the latter models, in particular #4a 
and #4b, to stellar mass loss. For these models mass loss is not adiabatic as 
is the case for simulations #1 and #2, but rather impulsive (see (T^J for a 
discussion in relation to variations in the orbital parameters of binary star 
clusters) . The later increase in the half-mass relaxation time is mainly caused 
by stellar mass loss and, in a lesser extend by the internal structural changes 
resulting from the internal dynamical evolution. The relaxation time for the 
various simulation models tend to evolve in almost parallel trajectories, with 
small variations. 

We approximate the time evolution of the relaxation time with the follow- 
ing two equations, for the N = 12k clusters we adopted 



Here k = 1/10 for the 12k models and k = 1/6 for the 128k models. The 
resulting tracks are presented in Fig.[5] 

Since the evolution of the half-mass relaxation time for each set of simu- 
lations with the same number of stars behaves quite similar, with an initial 
offset, we decided to invert Eq.Q]to enable an extrapolation of the initial re- 
laxation time for the observed clusters. The result for the 21 observed clusters 
with an age 300 Myr is presented in Fig.0 Here we plot our estimate for 
the initial half-mass relaxation time (t*, ) as a function of the measured age of 
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Fig. 1. Evolution of the two-body relaxation time at the half mass radius for the 
four clusters #1 (bottom), #2, #3 and #4 (top). The solid curves represent the data 
from models #{1..4}a, dashes give #{1..4}b and dotted curves gives the results for 
models #{1..4}c. Model #2a was extended to an age of 1 Gyr to demonstrate that 
the generally behavior doesn't suddenly change drastically. 

the observed cluster. We ignore the cluster mass in the reconstruction of the 
initial relaxation times for the observed clusters. From Fig.^ however, it may 
be clear that the effect of the number of stars in the cluster is substantial. We 
therefore opted for having the N — 12k models to provide an estimate of the 
upper limit and N = 128k for providing the lower limit to our estimate of the 
initial relaxation time. Those lower and upper limits are presented in Fig.E| 
as open and filled circles, respectively. 

4 Discussion and conclusions 

We performed extensive simulations of the young star clusters in the large 
Magellanic cloud. The simulations ignore the external tidal field, and we lim- 
ited the initial mass function to a Salpeter between 1M© and 100 M©. We 
realize that the adopted initial mass function may not be representative for 
the large Magellanic cloud, and we performed several additional simulations 
using the initial mass function proposed by J7j, which extends down to the 
helium burning limit. The difference with the results presented here, how- 
ever, are quite small, but wc tend to underestimate the initial relaxation time 
compared to using a more realistic initial mass function. 
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Fig. 2. Approximated evolution of the two-body relaxation time at the half mass 
radius for the four clusters #1 (bottom), #2, #3 and #4 (top). The solid curves 
represent the data from models with 12k stars (Eq.0 and dotted curves gives the 
results for models 128k stars. The bullets give an estimate for t T ix from the observed 
clusters by adopting a mean mass of 0.5 Mq and rh m = 3.92r CO re. 

The half-mass relaxation time for the simulated clusters evolves on almost 
parallel trajectories with an initial offset based on the initially selected re- 
laxation time. The trends in the relaxation time is somewhat different from 
the smaller (in TV) clusters than for the larger clusters, as is evident in Fig.^ 
We attribute this divergence to the difference in the the initial number of 
stars in these models, in particular since all other parameters were kept as 
much as possible identical. We note, however, that our simulation with 128k 
stars were computed with 10% primordial hard binaries, whereas the smaller 
simulations with 12k stars were computed without primordial binaries. This 
difference in the initial conditions may attribute to the differences in the evo- 
lution of the half-mass relaxation time. To test this hypothesis we performed 
several additional simulations with 64k stars and no initial binaries. Inter- 
estingly, the evolutionary tracks for the half-mass relaxation time for these 
models is roughly situated between the 12k and the 128k models. Based on 
these data we argue that the differences we observe between the 12k and the 
128k models can be attributed to the differences in the initial number of stars, 
and that the presence of primordial binaries has little effect. We will discuss 
these issues in more detail in an upcoming paper VI in the star cluster ecology 
series (Portegies Zwart, McMillan & Makino 2006+x {with x £l;i> 1}, in 
preparation). 
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Fig. 3. The initial half-mass relaxation time (tl lx ) as a function of the measured age 
of the observed LMC cluster. For each cluster we plot two symbols connected with a 
vertical line. The lower (open circle) symbol indicates extrapolation with time using 
Eq.0 for 128k whereas the upper symbols (bullet) gives the result from inverting 
Eq.0for 12k. 



Based on the here presented simulations we reconstruct the initial relax- 
ation time for young ( ^ 300 Myr) star clusters in the LMC with the data 
from . In Fig.0] we present the cumulative distribution of initial half-mass 
relaxation times for the observed star clusters in the LMC. For comparison we 
overlaid the distribution of the present day relaxation times for the observed 
clusters, which ranges over more than two orders of magnitude. We conclude 
that the initial relaxation time for the young ( ^ 300 Myr) star clusters in the 
large Magellanic cloud ranges from ~ 200 Myr to ~ 2 Gyr. 
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Fig. 4. Cumulative distribution of the initial half-mass relaxation time for the ob- 
served sample of star clusters in the LMC. The open circles represent the distribution 
when the initial half-mass relaxation time was reconstructed using the simulations 
with 128k stars, the filled circles are reconstructed using the 12k simulations. The 
solid curve without points (to the right) gives the observed distribution of present 
day relaxation times for the LMC clusters from data published by pQ. We only used 
the data for those clusters that are younger than 300 Myr. 
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